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Theoretical studies on charge ordering phenomena in quarter-filled molecular (organic) con- 
ductors are reviewed. Extended Hubbard models including not only the on-site but also the 
inter-site Coulomb repulsion are constructed in a straightforward way from the crystal struc- 
tures, which serve for individual study on each material as well as for their systematic under- 
standings. In general the inter-site Coulomb interaction stabilizes Wigner crystal-type charge 
ordered states, where the charge localizes in an arranged manner avoiding each other, and can 
drive the system insulating. The variety in the lattice structures, represented by anisotropic 
networks in not only the electron hopping but also in the inter-site Coulomb repulsion, brings 
about diverse problems in low-dimensional strongly correlated systems. Competitions and/or 
co-existences between the charge ordered state and other states are discussed, such as metal, 
superconductor, and the dimer-type Mott insulating state which is another typical insulating 
state in molecular conductors. Interplay with magnetism, e.g., antiferromagnetic state and spin 
gapped state for example due to the spin-Peierls transition, is considered as well. Distinct sit- 
uations are pointed out: influences of the coupling to the lattice degree of freedom and effects 
of geometrical frustration which exists in many molecular crystals. Some related topics, such 
as charge order in transition metal oxides and its role in new molecular conductors, are briefly 
remarked. 

KEYWORDS: molecular conductors, charge ordering, strongly correlated electron system, low-dimensional 
systems, metal-insulator transition, superconductivity, magnetism, electron-lattice coupling, 
geometrical frustration 



1. Introduction 

\ The charge ordering (CO) phenomenon is actively 

■ studied in the research field of charge transfer type 
[ molecular conductors,^ since it plays a key role in their 

physical properties. Following its clear observation in 
, DI-DCNQl2Ag,^ metal-insulator transitions in many of 

■ these materials are now understood as due to CO. It has 
' been found not only in newly synthesized compounds but 

■ also in systems which have been well known for years, 
\ where, however, its existence has been veiled until re- 
cently. Typical examples are two representative families 
of this field, the Bechgaard salts TMzX^ (TM = TMTSF 
or TMTTF) and ETsX^ (FT = BFDT-TTF) where X 
takes different atoms or molecules, both being studied for 

• more than 20 years. The CO transition was revealed re- 
cently and now became one of the most important issues 
in these systems. 

These materials are members of the so-called 2:1 salts 
expressed as A2B, to which interest of this field has been 
conducted. Numerous compounds with such a 2:1 com- 
position have been synthesized and found to exhibit a 
rich variety of properties. Many of them show elec- 
tron conduction at room temperature, where the car- 
riers are due to a charge transfer from cations S"*" or 
anions B~ , resulting in an average valence of —1/2 or 
-1-1/2 for A molecules, respectively. The B ion has closed 
shell in most cases then the valence band near the Fermi 
energy is composed of the frontier orbital, LUMO or 



HOMO, of the A molecule, which is quarter-filled as 
a whole in terms of electrons or holes. The variety in 
their properties has been revealed to be originated from 
the diversity of anisotropic lattices resulting in differ- 
ent non-interacting band structures, together with strong 
correlation effects experienced by electrons among this 
HOMO/LUMO band determining the low energy prop- 
erties.^ CO is a typical consequence of such strong cor- 
relation, namely, large electron-electron Coulomb repul- 
sion compared to the kinegic energy, especially due to 
the long-range nature of this Coulomb force. In fact, it is 
now ubiqitously found in A2B compounds as well as in 
other strongly correlated electron systems such as tran- 
sition metal oxides. 

In this article, we review theoretical aspects of CO 
mainly aiming at following points: In what kind of situ- 
ation are they formed? In what situation do they melt? 
How does the spin degree of freedom act? Are there any 
superconducting (SC) state near the CO phase? 

Such theoretical works on CO have been done from 
early days, motivated by experiments. For example, the 
metal-insulator transition in a classical transition metal 
oxide Fe3 04, the magnetite, was proposed to be due to 
CO by Verwey,^ although its existence is still contro- 
versial to date.^° Another trigger was an early molec- 
ular conductor TTF-TCNQ, where an incommensurate 
Ak-p charge-density-wave (CDW) is observed, but only 
in a diffusive manner therefore long ranged order is not 
achieved. Here we refer to the term CO as the phe- 



2 J. Phys. Soc. Jpn. 



Full Paper 



H. Seo, J. Merino, H. Yoshioka, and M. Ogata 



nomenon due to strong Coulomb interaction, sometimes 
called as "Wigner crystal on lattice". This should be 
distinguished with other transitions resulting in peri- 
odic modulations of the charge density, such as the 2fcF 
CDW (the Peierls-Frolich state) driven by the nesting 
of the Fermi surface together with the electron-lattice 
coupling, which is essentially a phenomenon at weak 
correlation. 

The observations of CO transition in A2B molecular 
systcims have stimulated many theoretical studies ad- 
hered to these compounds. Experiments have fortunately 
appeared around when researchers in this field started 
to realize that effects of electron correlation could be 
modeled in a straightforward way and that rather simple 
models would successfully describe their physical prop- 
erties.^^ That is, each constituent A molecule is repre- 
sented by a "site" and only the frontier orbital is con- 
sidered. The non-interacting band structures near the 
Fermi level are well reproduced by the extended Hiickel 
tight-binding scheme. ^^'^^ Then, the Coulomb interac- 
tion between electrons in this orbital is taken into ac- 
count. ^^■'^^ We call such model as the extended Hubbard 
model (EHM), written as follows: 

Wehm = - ^ (*u4o-Cjo- + h.c^ 

+ '^Uni^nii +'^Vijninj. (1) 

» {ij) 

Here, {ij) denotes pairs of the lattice sites (i.e., 
molecules) i and j, a is the spin index which takes t 
and I, riia- and cj^ (cjo-) denote the number operator and 
the creation (annihilation) operator for the electron of 
spin (T at the ith site, respectively, and rij = nj| + n^j. 
The transfer integrals, tij, reflect the anisotropy resulting 
from the particular spatial extent of the frontier orbital, 
calculated, e.g., by the extended Hiickel method or from 
tight-binding fitting of first principle calculations. The 
Coulomb interactions of not only on-site U but also inter- 
site Vij are considered, the latter being crucial for the CO 
as we will see later explicitly. In the [^2]~-B''' {[A2]'^B~) 
systems the non-interacting band as a whole is quarter- 
filled in terms of electrons (holes), namely, there exists 
one electron (hole) per two sites on average. 

Because of this clear way of constructing microscopic 
models from crystal structures, results of theoretical 
works could be checked back in the experiments, and such 
interplay has greatly accelerated the research. Through- 
out these we have learned that, although the basic pic- 
ture of CO is rather classical and essentially known from 
the early days, the physics therein is rich and diverse. 
The main reason for such diversity is the variety in the 
geometry of lattice structures of the materials where it 
is realized. More specifically, relative positions between 
molecules not only refiect directly on the anisotropy in 
tij controlling the band structure, but also affect drasti- 
cally the nature of the CO state itself through V^j . This 
is in contrast with the case of the Mott insulating state in 
half-filled systems where the driving force is the on-site 
term U, which is a character of the atom/molecule itself; 
in the case of the dimer-type Mott insulator in A2B sys- 



tems (see later), it is the "on-dimer" Coulomb repulsion, 

TT 17,19,20 
*^dimer* 

The research of CO is still continuing and rather 
rapidly growing. New phenomena are uncovered, even in 
the above mentioned compounds, and now under exten- 
sive investigations. For example a pressure-temperature 
phase diagram of DI-DCNQl2Ag has been explored 
where an anomalous temperature dependence of metal- 
lic resistivity p oc is found just beyond the border of 
the CO phase. A peculiar interplay between CO and 
magnetic properties found in TMTTF2X compounds un- 
der pressure^^ requires a reformulation of the generic 
phase diagram of TM2X.^^ Problems of CO system on 
anisotropic triangular lattice structures characteristic of 
ET2X compounds have many new aspects.*' We can- 
not offer comprehensive explanation for each of these 
cases here as many are still not yet resolved theoreti- 
cally, but we hope that this review, by explaining the 
present status, would provide a base for tackling such 
new physics and lead the readers toward challenges. 

The organization of this paper is as follows. Since most 
of the A2B compounds have low-dimensional structures, 
we will devide this review into one- dimensional (ID) and 
two-dimensional (2D) problems, while in reality there 
exist finite interchain/interlayer interactions; they are 
quasi-lD and quasi-2D systems. Theoretical results de- 
voted for the quasi- ID compounds are described in § 2. 
This starts with studies on purely ID electronic models 
and then additional effects, such as interchain interac- 
tion and coupling to the lattice degree of freedom, are 
considered. The ID case can be studied in more con- 
trolled ways, analytically and numerically, than in the 
study of the quasi-2D compounds, which we discuss in 
§ 3. There, theoretical works on 2D electronic models 
are still in progress and influence of additional effects, 
e.g., coupling to the lattice, is not fully understood yet. 
Studies aimed at SC states near the CO phase in 2D 
models, motivated by its observations, are reviewed as 
well. A problem of CO systems under geometrical frus- 
tration which is expected to be relevant to many molec- 
ular conductors will be pointed out in § 2 and 3. In § 4 
related topics will be mentioned, such as analogous CO 
states observed in transition metal oxides. Possible roles 
of CO in other molecular systems will be added as well, 
as perspectives. A summary is given in § 5. 

Some of the experimental studies on CO are mentioned 
in this review but many references including important 
ones are left out; we refer to other reviews from exper- 
imental standpoints"'^' which would be complementary 
to this article. One can find many review articles on the 
properties of molecular compounds in general,^' ^"^ espe- 
cially we refer to refs. 8 and 18 for papers from theoreti- 
cal but more systematic point of views including the CO 
systems as well. 

2. Quasi-One-dimensional Systems 

Early theoretical works on CO in ID models have been 
performed motivated by the observation of Akp CDW 

in TTF-TCNQ as mentioned in § 1, where the impor- 
tance of the long-range Coulomb interaction was empha- 
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sized. However for the quasi-lD A2B systems, be- 
fore CO was found, analyses were mainly concentrated 
on the Hubbard-type models only considering the on-site 
Coulomb energy U, but some did discuss the relevance 
of the inter-site Vij.^^^^^ For example, Mila^'^ estimated 
that the nearest-neighbor Coulomb repulsion is appre- 
ciable in TM2X, such as more than one third oi U, by 
comparing calculations on the dimerized version of ID 
EHM (see eq. (2)) with optical measurements.'^'* After 
CO was found, many further works devoted to the ID 
EHM and its variants have been carried out, which we 
will review in this section. 

2.1 Charge ordering in quasi- one- dimensional systems 

When the first observation of CO in DI-DCNQIzAg 
was made by ^^C-NMR,^ in an independent theoreti- 
cal work based on mean-field (MF) approximation, Seo 
and Fukuyama'^^ proposed that CO due to the nearest- 
neighbor Coulomb repulsion might exist in TMTTF2v'C. 
This was based on a comparison between self-consistent 
solutions at zero temperature obtained by the standard 
MF treatment applied to the appropriate dimerized ID 
EHM and the spin structure in the antiferromagnetic 
phase suggested by a -^H-NMR measurement."^^ Soon af- 
ter, it was found that members of TMTTF2A" indeed 
show CO, directly seen in a ^^C-NMR measurement.'^ 
Moreover a divergence of the dielectric constant at the 
CO transition temperature is observed, suggesting a fer- 
roelectric state. "^"^ This is a consequence of CO, which we 
will mention later in this section. 

These systems are members of DCNQl2Ar {X: mono- 
valent metal cation X+, e.g., Ag and Li)^^'^^ and TM2X 
{X: monovalent anion X~, e.g., PFe, AsFe, SCN, and 
23 respectively, both having quasi-lD structures. DC- 
NQI stands for the i?ii?2-DCNQI molecule where Ri,R2 
are sustituents such as CH3, Br, I, etc. Here, quasi- 
lD "structure" implies twofold meanings: in their crys- 
tal structures the DCNQI/TM molecules assemble in a 
stacking manner, and, in their electronic structures the 
transfer integrals in the interchain direction, Winter 1 sn'e 
one order of magnitude smaller than those in the intra- 
chain direction, iintra- 

The networks of DCNQI/TM molecules are schemati- 
cally shown in Fig. 1. In DCNQl2Ar the interchain cou- 
plings are three-dimensional (3D) and uniform, while 
TM2Ar they are quasi-2D and the TM molecules are con- 
nected in a rather complicated way (the /?-type structure; 
see § 3.1). For DCNQI2X systems, |ii„tra| ^ 0.15 ~ 0.25 
eV, Itintorl ^ 0.01 ~ 0.03 eV,''°''^i while for TM2X sys- 
tems: in TMTSF compounds |ii„tra| ^ 0.2 - 0.4 eV, 
Ifinterl ^ 0.01 - 0.05 cV and in TMTTF compounds 

Itintral - 0.1 - 0.25 eV, liinterl ^ 0.01 - 0.03 eV.lf^'42-44 

These values depend on the actual salts and different 
methods of calculation also provide varied estimations. 

We note that the distances between molecules does 
not simply correspond to the degree of anisotropy in 
tij, since the anisotropic shape of the frontier orbitals 
makes the dependence of tij to the relative configuration 
of molecules rather complicated in general. *^'^^ In con- 
trast, Vij obeys more or less a monotonic function of the 




Fig. 1. (Color online) Schematic representation of the molecular 
networks in (a) DCNQI2X and (b) TM2X. The crystallographic 
axes are shown; in (a) we take a = a + b,b = a — b where a and 
b are the axes in the tetragonal cell. The interchain interactions 
are three-dimensional in (a), whereas it is two-dimensional in (b) 
since interlayer couplings are one order of magnitude smaller due 
to the anion layers. The indices (for TM2X, taken from ref. 16) 
are not only for tij but also for different values of Vij , although 
the degree of anisotropy can be different (see text). 

distance because it is the Coulomb interaction. There- 
fore the interchain part of Vij is not necessarily small 
compared to the intrachain one even in these quasi-lD 
systems.*^' However, let us neglect the interchain in- 
teractions first and discuss this issue later in § 2.4 and 
2.6. 

Estimations of the Coulomb energies in these molecu- 
lar systems are difficult at present. Quantum chemistry 
calculations are performed for an isolated molecule or 
clusters of them, which provide unrealisticly large values 
since the screening effect in solids is left out, while such 
estimations from the first principle in these molecular 
solids are still yet to be done. However it is believed to 
be of the order of [/ ~ 1 eV from different measurements 
and estimates, which gives J7/|iintra| — 5 for DCNQl2Ar 
and TMTTF2X, while for TMTSF2A: smaller values of 
C^/I^intral — 3. The estimated values for the intrachain 
Coulomb energy T^ntra are again ambiguous but many 
provide rather large values: the ratio l^ntra/C^ in a range 
about 0.2 - 0.6.33' 46-48 

A crucial difference between these two families is that 
the stacking of DCNQI molecules is uniform while that 
of TM moluecules is slightly dimerized, as seen in Fig. 1. 
Thus a minimal effective model to investigate the occu- 
rance of CO in these compounds is the ID dimerized 
EHM, represented as, 

Um = -t^(l + (-l)'<5d)(cl+i,c,. -f h.c.) 

i(T 

+U^ni'^nii +V^nini+i, (2) 

i i 

where i is the site index along the chain. Here, the trans- 
fer integrals allow dimerization as alternating t{l + S^) 
and t(l — (5d); in the DCNQI compounds 5d = 0, while for 
the TM compounds values of tai and ia2"'^^''*^"^'* (for the 
indices see Fig. 1) read <Q.l. The inter-site Coulomb 
repulsion between neighboring sites is set to be uniform 
as V , which is an approximation for TM2Ar. Again, the 
value of 6(1 does not directly result in a similar value of 
dimerization in Vintra- In fact, quantum chemistry calcu- 
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Fig. 2. (Color online) Two limiting cases of strongly correlated 
insulators in quarter-filled system*' : (a) the charge ordered 
insulating state and (b) the dimer Mott insulating state. Black 
dots and colored area represent the lattice sites and the localized 
carriers, respectively, while the thickness of the bonds show the 
difference in the transfer integrals. The arrows show the param- 
agnetic localized spins. 

lations for TM2X''^'"'^ provide the dimerization parame- 
ters in Vai and Va2 to be less than 1 %. For the DCNQI 
compounds the electron filling is one quarter while it is 
three quarter for the TM compounds, but in this ID 
model the two situations are equivalent since electron- 
hole symmetry holds. Note that other effects may break 
this symmetry then we should treat the two cases sepa- 
rately. 

It is useful to describe two limiting cases for the insu- 
lating states in this model at quarter-filling due to strong 
Coulomb interaction,^'^* which are shown in Fig. 2. 
These states, in a broad sense, are realized in a wide 
range of A2B compounds, not only in quasi-lD but also 
in quasi-2D systems discussed in § 3. One is the Wigner 
crystal-type CO state stabilized by V in the presence 
of strong U, and the other is the Mott insulating state 
stabilized by U in the presence of strong dimerization 
6d' a dimer Mott insulator. The spin degree of freedom 
in such insulating states would behave as an = 1/2 
localized spin chain, where the spins are located on ev- 
ery other site for the CO state whereas on every dimer 
for the dimer Mott insulating state. Roughly speaking, 
charge localization is determined by large energy scales 
such as U and V, while spin properties are determined 
by a smaller energy scale of the order of the Heisenberg 
coupling, J, acting between these localized spins. We will 
discuss their detailed properties in the following subsec- 
tions. 

The MF solutions mentioned above are consistent with 
such two limiting cases. This is seen in the obtained 
spin and charge patterns for the two cases, schematically 
shown in Fig. 3, and actually many works have been 
performed on Hid based on such MF results. How- 
ever we should keep in mind that this MF treatment 
cannot correctly describe such insulators at strongly cor- 
related regime in general, for example the paramagnetic 
insulating phase at temperatures above the magnetically 
ordered phases observed in experiments cannot be repro- 
duced. Moreover in purely ID models the role of quantum 
fluctuations which is left out in MF is crucial, and we will 
see that it considerably modifies the MF results. 

We note that the cations/anions are at crystallo- 
graphically equivalent positions from the DCNQI/TM 
molecules at high temperatures so that, essentially, they 
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Fig. 3. (Color online) Typical mean-field solutions in one- 
dimensional dimerized extended Hubbard model at quarter- 
filling,'^^ WiD in eq. (2): (a) the charge ordered antiferromag- 
netic insulating state for <5(i = and large U and V, and, (b) the 
dimer-type antiferromagnetic insulating state for 5^ 7^ 0, large 
U, and V = 0. The size of the colored circles and the arrows 
represents the charge density and the amount of spin moment 
on each site. The charge densities are (a) disproportionatcd al- 
ternatively as, 1/2-1-5, 1/2-5, 1/2-1-5, 1/2-5, where 5 is the 
amount of charge disproportionation, and (b) uniform at 1/2. 

do not contribute to the electronic properties. However, 
in the case of non-centrosymmetrical anions in TM2X, 
the anions show a disorder-to-order transition by lower- 
ing temperature and generate a potential with longer pe- 
riodicity than the original unit cell, resulting in modifica- 
tions in the one-particle properties. ^° We do not discuss 
such cases in this paper since its relation with the CO 
due to electron correlation is still obscure, but the fact 
that the anion ordering produces charge density mod- 
ulations suggests the importance of the anions in some 
cases. Actually possible roles of the X unit in stabilizing 
CO states through its coupling to the electron system 
will be addressed in § 2.5. 

2.2 One- dimensional extended Hubbard model 

The quarter-filled ID EHM, Hid in eq. (2) with 6d=0, 
was first considered by Ovchinnikov^^ who discussed that 
in the U/t = 00 limit a CO insulating ground state is re- 
alized for V > 2t {= Vcr), as in the following. In this limit, 
double occupancy is strictly prohibited so that every site 
is either occupied or unoccupied. Furthermore the spin 
degree of freedom is frozen since the "order" of electrons 
cannot be changed due to the ID nature. Then the charge 
degree of freedom is equivalent to a spinless fermion (SF) 
model at half-filling with the nearest neighbor Coulomb 
interaction V, i.e., the interacting SF model, and the spin 
degree of freedom acts freely as Curie spins. 

This SF model, by identifying occupied site as up spin 
and imoccupied site as down spin, can be mapped onto a 
ID 5 = 1/2 (pseudo-)spin model through Jordan- Wigner 
transformation. Quarter-filling in the original electronic 
system corresponds to half-filling in the SF model where 
half of the sites would be occupied, and to total magnetic 
moment of zero in the effective spin model. The inter-site 
Coulomb interaction between SF transforms to an anti- 
ferromagnetic interaction, J^SfS^ with — V, while 
the kinetic enery term becomes as J^y^Sf SJ+Sf Sj) with 
Jxy = 2t. This is an XXZ model, where a T — Q phase 
transition from a gapless "XY" state (J^ < J^y) to an 
antiferromagnetic "Ising" state (J^ > Jxy) is known. ^-'^ 
Transforming back to the SF model, these correspond 
to a metallic Tomonaga-Luttinger liquid (TLL) and the 
CO insulating state, repectively; the critical point is the 
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Heisenberg point Jxy = Jz, i-e-, V = Kr- 

It is noteworthy that this is in fact a rigorous demon- 
stration of a "classical" picture of CO discussed in sys- 
tems regardless of the dimension. In general, the CO 
problem, in the condition of kinetic energy set to 
{tij = 0) in addition to infinite U with no double oc- 
cupancy discussed above, can be described by SF's in- 
teracting only by the inter-site Coulomb repulsion Vij. 
This can be mapped onto Ising models, JijaiCfj, with 
Jij = Vij, similarly to the J^-term in the ID case above. 
Therefore, analogies between spin systems and CO sys- 
tems can be generally expected in the presence of strong 
correlation, although the mapping is only approximate 
for 2D and 3D systems with finite tij and U, and even 
for ID systems with finite U. 

Since the charge and spin degrees of freedoms are com- 
pletely decoupled in the U/t — co limit, the ground state 
wave function has the form |$) = |</>sF)|Xcr), where |(^sf) 
is the ground state of SF's and \xa) denotes the spin 
configuration. On the other hand, at finite U/t, these 
spins interact with each other. In the large-?/ region, 
recently Tanaka and Ogata®^ analytically discussed the 
F-dependence of the magnetic susceptibility x- This is 
done by extending the \axge-U studies of the ID Hub- 
bard model based on the Bethe ansatz,^''' ''''' where the 
charge degree of freedom can be represented by SF's, 
and the spin degree of freedom is described by Heisen- 
berg spin chains with an effective exchange coupling Jeff. 
For the ID EHM, according to the degenerate perturba- 
tion theory, the degeneracy of the spin system is lifted 
by a term proportional to J = At^ /U. However, because 
of the presence of the charge degree of freedom, J^s is 
not just equal to J, but also depends on the probabil- 
ity of finding two SF's on nearest-neighbor sites, namely, 
(nini+i)sF = (0SF|?ii"i+i|?!'SF)- (There is another con- 
tribution to Jeff discussed in ref. 53.) This value de- 
creases as V increases and its derivative with respect 
to V is continuous even at Vcr/t,^^ where the CO tran- 
sition occurs. This is because the phase transition is a 
Berezinskii-Kosterlitz-Thouless (BKT) type,'"'^' '''^ which 
is essentially different from either first or second order 
transition. The divergence of the correlation length is 
more rapid than any power law, and then the gap in- 
creases very slowly near the critical point: the charge 
gap behaves as Ac ~ 47rexp(-7r^/2(2(p-l))^),^^'^^ with 
p = V/2f. > 1. Reflecting these, x increases as a mono- 
tonic function of V since x is proportional to l/Jeg,^'*'^^ 
continiously through the CO transition point (see also 
Fig. 5). This indicates that the fluctuation is large near 
the phase boundary even in the CO states, and conse- 
quently CO has small effect on the spin degree of free- 
dom. 

For whole finite values of U and V, Mila and Zotos^^ 
first provided the ground state phase diagram on the 
plane of U/t and V/t, using numerical Lanczos exact di- 
agonalization (ED) method^^ for finite size systems up to 
L = 16 sites. They have judged the metal- insulator phase 
boundary by calculating the charge gap for finite sys- 
tems as Ac = [Eo{L/2 + l) + Eo{L/2-l)-2Eo{L/2)]/2, 
where Eo{ri) is the ground state energy for n electron 
system (n = L/2 for quarter-filling), and extrapolating 




U/t 

Fig. 4. Ground state phase diagram of the one-dimensional ex- 
tended Hubbard model at quarter-filling on the plane of U/t 
and y/t.®^ CO stands for charge order and TLL for Tomonaga- 
Luttingcr liquid. SC represents the superconducting phase. Con- 
tour map for the Tomonaga-Luttinger parameter Kp is shown. 
The bold line represents the boundary for the metal-insulator 
transition. The shaded area is the region with an exponentially 
small gap. [By courtesy of S. Ejima.] 
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Fig. 5. Spin susceptibility x calculated by exact diagonalization 
as a function of Y for several values of U at quarter filling (t = 
1).^^ Solid squares represent exact results in the one-dimensional 
Hubbard model and arrows indicate the points above which the 
charge gap opens. [After ref. 53.] 

the results to the thermodynamic limit. Later, differ- 
ent numerical methods have been applied by different 
authors'^^'^^"^'^ and the main features are confirmed: 
the phase diagram by Ejima et al^'^ is shown in Fig. 4, 
which is determined by the density matrix renormaliza- 
tion group (DMRG) method treating system sizes up to 
more than L — 100 sites. The CO insulating state is 
in the large {U jt, V/t) region, while the metallic state 
in the rest is the TLL. The TLL parameter Kp, which 
is the exponent characterizing the power of the correla- 
tion function, becomes equal to Kp = 1/4 at the phase 
boundary {Kp = 1 at [/ = V = 0). As /7 is increased, 
this phase boundary approaches to [U/t = 00, V/t = 2), 
the exact value explained above. 

The magnetic susceptibility x has been discussed in 
ref. 53 for such finite {U/t, V/t). This is done numeri- 
cally by the Lanczos ED method up to L = 16 sites, as 
shown in Fig. 5 where x at T = for several values of 
U/t and V/t are plotted. The TLL property x = 2/'ku„ 
is used with being the velocity of the spin excita- 
tion (i.e., the spin velocity). The critical values for the 
metal-CO insulator transition are indicated by arrows in 
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Fig. 6. (Color online) Coexistence of charge order and dimeriza- 
tion in the one-dimensional dimerized extended Hubbard model 
for TM2X. This is a schematic drawing of the mean-field solu- 
tion, where antiferromagnetism also coexists.'^® Charge dispro- 
portionation exists as in the charge ordered state in Fig. 3 (a) 
whereas the amount of spin moment, parallel whithin dimers 
and antiparallel between them, are also disproportionated in a 
similar way. 

Fig. 5. One can see that x increases with V , namely, Jeff 
is suppressed as in the large-J7 case above, x smoothly 
varies as a function of V even through the CO phase tran- 
sition, as in the exact calculation in the large-?/ case.^'^ 
The spin susceptibility is again not affected much by the 
opening of a charge gap. This is consistent with field the- 
oretical studies that we will see in § 2.3, indicating the 
transition to be of the BKT type for finite (C//i, V/t) as 
well. 

In the presence of finite dimerization (5d, it has been 
well known for the V — Q case, namely, the ID dimer- 
ized Hubbard model, that together with U always give 
rise to the dimer Mott insulating state (Fig. 2(b)) for 
quarter-filling. This is because 5^ opens an energy gap 
at ±2A;f in the non-interacting energy dispersion making 
the lower band effectively half-filled. Then even for an in- 
finitesimal repulsion the system becomes a Mott insula- 
tor which is characteristic of ID half-filled electronic sys- 
tems. This situation have been theoretically investigated 
extensively,^^ from the large-C/ case^^'^^ to the small-t/ 
case.^* When finite V is added, there arises a competition 
between two insulating states, this dimer Mott insulator 
and the CO insulator. The CO state here can be con- 
sidered as a coexistence of dimerization and CO as seen 
in Fig. 6. Accurate numerical calculations are yet to be 
done and now in progress, ^'^'^^"^^ so let us discuss such 
competition in the next subsection. 

2.3 Bosonization 

The bosonization method is one of the most powerful 
tools to treat ID systems analytically.^^ In this approach, 
quantum fluctuation, which plays an essential role on the 
electronic state in ID systems, can be fully taken into 
account, in contrast to the MF theory. Furthermore, the 
analytical treatment enables us to capture insight of the 
physics, e.g., found in numerical calculations discussed 
in the previous subsection. It has been successfully ap- 
pHed to Hid in eq. (2) by Yoshioka et al.J^''^'^ as in the 
following, who clarified the relation between such lattice 
models and the conventional "g-ology" picture. 

In the ordinary bosonization procedure, only the one- 
particle fermion states around the two Fermi wavenum- 
ber, ztfep, are considered and the four body interaction in 
the Hamiltonian is expressed in terms of these states. In 
quarter-filled systems, the Sfcp-Umklapp scattering ex- 
ists ''^"^^ because of kp = n/{Aa) with a being the lat- 
tice spacing and it is crucial for the appearance of the 
insulating state. However, one can not obtain this Umk- 
lapp scattering by the ordinary method above. This is 
because it can only be expressed by the interaction pro- 
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Fig. 7. (Color online) Diagrams representing the Sfcp-Umklapp 
scattering. Here p = +/— expresses the right/left-going one- 
particle states and the solid and the dotted lines express the 
electrons near ztfcp and ±3fcF, respectively. 

cesses where four right-going electrons are scattered into 
the left-going states and vice versa to gain (or loose) Sfcp , 
and then the scattering should include higher order in- 
teraction processes through the one-particle state around 
±3/cF- 

A systematic way to derive the Sfcp-Umkpapp scatter- 
ing has been developped^^' as follows. The one-particle 
states are divided into two parts: the states near ifep and 
those near ztS/cp- The effective Hamiltonian written in 
terms of the former states near ztfep is obtained by inte- 
grating out the one-particle states near ztSfep and treat- 
ing the interaction processes including both the states 
near ±kp and those near ±3/cf perturbatively. The Skp- 
Umklapp scattering appear in the third order interaction 
processes shown in Fig. 7. 

Then we can express the effective phase Hamiltonian 
for the ID EHM with = 0, by using the bosonic phase 
variables representing the charge and spin fluctuations.^^ 
The Hamiltonian is divided into the charge part. Tip, 
and the spin part, Ti.^, which is essentially the same as 
that of the ID isotropic Heisenberg model with a gapless 
excitation. Tip is given by the sine-Gordon model shown 
below. 
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dx cos 4:9 „ 
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where [9p{x),—dx'(pp{x')/{2Tr)] = iS{x — x'), and Vp and 
Kp are the velocity of the charge excitation and the TLL 
parameter introduced in § 2.2, respectively. The quantity 
is the ultraviolet cutoff of the order of a~^. 
The coefficient of the non-linear term gi/4 express- 
ing the Sfcp-Umklapp scattering is written as, gi/4 = 
{Uaf{Ua - 4Fa)/8(7rto)2. The order parameter of CO 
at the i-ih site, Oco{xi) is given by using the phase vari- 
ables, as Oco{xi) oc cos(4fcpa:i +29p) = (— I)'cos20p, 
where Xi ~ ia. This indicates that CO is stabilized 
for 9p — and 9p ~ 7r/2, whereas it disappears for 
9p = 7r/4 and 9p — 3tt/A. The former and the latter 
can be realized in the case of f/1/4 < and gi/4 > 
0, respectively. This shows explicitly the fact that the 
inter-site Coulomb repulsion V stabilizes the CO state 
through the decrease of gi/4 < as y is increased. 
We note that this condition < for the appear- 
ance of CO, i.e., 4V > U, is also a necessary condition 
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in the MF treatment. This is based on the charge sus- 
ceptibiHty within random phase approximation (RPA) , 
Xciq) ^ Xoiq)/[i + (U + 4V^cosga)xo(g)], where xoiq) is 
the noninteracting charge susceptibihty taking a positive 
value; 4V > U is required for Xc to diverge at qa = w, 
i.e., for the realization of the CO state with q = Akp- 

The low energy property of the phase Hamiltonian Hp 
is determined by the renormalization group equations by 
introducing I — ln(a'/a) with a new scale a' (> a), 

^^Kp{l) = -^G„SfKp{l)\ (4) 

^Gi/4(/) = (2-8i^,(0)Gi/4(0- (5) 

The initial conditions are -K'p(O) = Kp and 01/4(0) = 
§1/4^1 {2iTVp). Whether the system is in the metalhc TLL 
state or in the CO insulating state is determined by 
the / — » 00 solutions characterized by Gi/4(oo) = or 
Gi/4(oo) = —00, respectively. The phase boundary de- 
termined by the renormalization group equations, where 
Kp{oo) = 1/4, is shown by the dotted curve in Fig. 8 
where a = laj-n is used. This phase transition is of BKT 
type since the sine-Gordon theory is equivalent to the 
low-energy properties of the classical 2D XY model. The 
phase diagram obtained by the present bosonization the- 
ory is qualitatively the same as that shown in Fig. 4. 

The bosonization theory above can be straightfor- 
wardly extended to include the dimerization together 
with finite V ^"^ As a result of the energy gap at ±2A;f in 
the non-interacting energy dispersion due to finite value 
of ^d, in addition to Hp in eq. (3), the 4fcF-Umklapp scat- 
tering is generated.^® This is expressed in terms of the 
phase variables as, 

ny2 = -,S^ [ dxsm20p, (6) 

where 171/2 is proportional to for \S(\\ <C 1. The ground 
state phase diagram in the presence of determined by 
the renormalization group equations, on the plane of V/t 
and U/t, is shown in Fig. 8. Due to the 4fcF-Umklapp 
scattering, the metallic TLL state changes to the dimer 
Mott insulating state with uniform charge distribution 
even for infinitesimal interactions, as expected from the 
works for V = 0.^^ The CO insulating state (Fig. 6) is 
suppressed because the commensurability energies due to 
the two kinds of Umklapp scattering compete with each 
other for 51/4 < and 7ii/2 has a larger scaling dimen- 
sion. Nevertheless, importantly, the CO phase still exists 
in the large {U/t, V/t) region. This model for the charge 
sector having two non-linear terms is called the double 
frequency sine-Gordon model and it is suggested that the 
phase boundary between the dimer Mott insulator and 
the CO insulator is an Ising type phase transition, where 
the charge gap becomes at the phase boundary.''^' 

2.4 Interchain coupling 

In purely ID systems, phase transitions do not occur at 
finite temperatures because any ordered state collapses 
due to thermal fluctuations. In this sense, the critical 
temperature for CO, Tco in the ID EHM that we have 
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Fig. 8. Ground state phase diagram of the one-dimensional 

dimcrizcd extended Hubbard model at quarter-filling on the 
plane of V/t and U/t, determined by the renormalization group 
equations for the bosonized phase Hamiltonian. The solid curve 
expresses the phase boundary between the charge ordered (CO) 
insulator and the dimer-typc Mott insulator for = 0.1,**''' while 
the dotted curve shows the boundary for = 0.0 where the Mott 
insulating phase is replaced by the Tomonaga-Luttinger liquid^^ 
(see Fig. 4). 

discussed in § 2.2 and 2.3 is always absolute zero temper- 
ature. However, quasi-lD materials such as DCNQI2X 
and TM2X undergo phase transitions at finite tempera- 
tures, which is a consequence of the three (or two) dimen- 
sionality. In the electronic sector, there exist two kinds 
of interchain couplings; one is the single particle hopping 
between the chains (tintcr) and the other is the mutual 
Coulomb interaction between electrons in different chains 
(Vinter)- As mentioned in § 2.1, fintra/iinter is of the order 
of 10 in the DCNQI and TM salts, while the interchain 
Vinter uiay bc of the same order with the intrachain Vlntra- 
There are quantum chemistry studies for TM2X actually 
showing that it the case."*^'^^ 

The interchain hopping Winter generates warping of the 
Fermi surface in the non-interacting band structure. In 
TMTSF2X compounds, because of this effect, together 
with the weaker effect of Coulomb repulsion than in their 
TMTTF analogs, the strongly correlated insulators such 
as CO and dimer Mott insulators arc not realized. Then 
the Fermi surface persists down to low temperatures. 
Such a dimensional crossover in half-filled models are the- 
oretically discussed^^ while that in the quarter-filled CO 
systems remains to be investigated. This should be rele- 
vant to the pressure- induced melting of CO observed in 
DI-DCNQl2Ag,^-'^ where the CO transition turns from a 
second order to a first order phase transtition by applying 
pressure, and a region with T^-resistivity appears in the 
metallic state at the verge of the CO phase. In TM2X, in 
contrast, the system becomes unstable due to the nesting 
of the Fermi surface toward a metal-insulator transition 
to an incommcnstirate spin-density-wave (SDW) ground 
state, which have also been studied intensively.^' -^^'^^ 

Effects of the interchain Coulomb repulsion on the CO 
state have been discussed recently by Yoshioka et al.^^ by 
using the bosonization method described in the previous 
subsection. When one considers only one kind of uniform 
coupling in the interchain direction, it is expressed as, 

W_L = Vj. ^ ni,;ni,i+i , (7) 
where I denotes the chain index. This may be applied 
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to the DCNQI compounds as the interchain couphng is 
uniform and isotropic while in TM2X the network is 
more compUcated as seen in Fig. 1 (we will consider 
them in § 2.6). In the "interchain MF" treatment,^'' 
the interchain part is treated within MF approximation 
as, Ui^irii^i+i {ni^i)nt^i+i + ni^i{ni^i+i) - {ui^i) {riij+i) , 
while quantum fluctuation within the ID chains is in- 
corporated. When (rij) along the chains is written as 
1/2 + {—lyS where S being the amplitude of CO, and an 
antiphase pattern of the Wigner crystal-type CO is as- 
sumed, the problem is reduced to an effective ID model 
where 5 is determined self-consistently. If S is finite, an 
energy gap opens at ±2A;f and as a result the ikp- 
Umklapp scattering, Hi/2) is generated, similarly to the 
case with finite dimerization discussed in § 2.3. However 
this Umklapp scattering is written with a cosine function 
of the phase variable but not with a sine function as in 
eq. (6), as follows, 

2(7ra)^ J 

where gi/2 is proportinal to S as long as |5| <C 1. There- 
fore, different from the case with dimerization the gener- 
ated non-linear term does not compete but cooperate 
with the cos46'p term when gi/4 < in eq. (3). We 
note that such term, which acts as an alternating on-site 
potential upon the ID chain, has been discussed from 
another context. Thus, the interchain interaction sta- 
bilizes the antiphase CO, and the critical temperature 
Tco becomes finite. Renormalization group study on this 
phase Hamiltonian*'^ indicates that an infinitesimal value 
of V± gives rise to finite Tco as long as Kp < 1/2 in 
the original ID EHM model. This results in a sudden 
enlargement of the CO insulating region in the U — V 
phase diagram once V± is considered, from the case of 
the ID EHM where the phase boundary is determined 
by Kp = 1/4, as in Fig. 4. 

2.5 Coupling to the lattice 

ID electron systems arc unstable when there exist cou- 
plings to the lattice degree of freedom. A well known case 
is the 2fcF CDW (Peierls) instability, where electron gas 
becomes unstable at low temperatures toward a lattice 
distortion of period 2A;f when the nesting condition of 
the Fermi surface is satisfied. In the purely ID case 
this condition is perfect so that any ID electron system 
is susceptible to the Peierls-Frolich state. Another case 
is when the electrons are localized due to strong correla- 
tion and the spin degree of freedom is described by a ID 
Heisenberg chain. There, the spin-lattice coupling makes 
the system unstable toward a spin-singlet formation ac- 
companied with a lattice distortion, i.e., the spin-Peierls 
(SP) transition. ^^'^^ In fact, SP transition is frequently 
observed in many quasi-lD A2B compounds, while its 
interplay with the charge sector such as the CO and the 
dimer Mott insulators is not obvious. 

Let us first consider the so-called Peierls-type coupling 
to the lattice degree of freedom, treated as classical vari- 
ables here neglecting quantum fluctuations, in addition 
to HiB in eq. (2) with 6d = 0. This allows modulations 
of the transfer integrals at the cost of the elastic energy, 
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i.e., t — > t{l + Ui^i+i) + Kuf ^_^_l/2, where Ui^i+i is the di- 
mensionless modulation of the lattice, controlled by the 
actual movement of the molecules i and i + 1. K is the 
renormalized lattice constant, which is taken to be uni- 
form here but in general, as in the discussion in § 2.1 on 
Vij, can be modulated; for example in the ID model of 
TM2X it would be dimerized as Kai and Ka2 (sec Fig. 1). 
Hereafter we write = Ui for simplicity. The "soft- 

ness" of the lattice is indicated by l/K; at ^ 00 the 
system is in the "hard" limit so that the lattice does not 
move and then the model reduces to the purely electronic 
model. Hid- 

One serious effect of this coupling is that it can gen- 
erate the dimer Mott insulating state even without in- 
trinsic dimerization {Sd = 0). This is understood in the 
U/t = 00 limit discussed by Bernasconi et al.^^ Extend- 
ing the discussion in § 2.2, the charge degree of freedom 
described by the half-filled interacting SF or equivalently 
the S* — 1/2 XXZ model is now coupled to the lattice. In 
the latter model the Peierls coupling appears in the XY 
term, J^y] for J^y = 2t > — V , infinitesimal spin lat- 
tice coupling l/K can drive the ground state toward SP 
lattice dimerization.^^ This is actually the dimer Mott 
insulator when one transforms back to the original elec- 
tron model. This state competes with the CO state at 
large V/t, and the transition point Vcr/t increases with 
increasing l/K.^"^'^^ 

When U/t is finite, the spin degree of freedom 
(in the original electronic model) becomes active. The 
model have been treated by different numerical tech- 
niques^^' ^^"^^ as well as by the bosonization method. ®^ 
As discussed in the previous subsections, the latter treat- 
ment can clarify the physical origin of each state which 
is realized. A phase Hamiltonian is derived following the 
treatment in § 2.3, where nonlinear terms additional to 
Hp in eq. (3) and Ha are generated in the presence of 
the lattice dimerization Ud and tetramization Ui, with a 
cost of lattice elastic energy Hei- These can be expressed 
as 

Hd = -^ [dxsin29p, (9) 
2(770;)^ J 

= "■:T7^~T2 / da;sin(^p - xt) cos9a, (10) 
2[Tray J 

when the lattice distortion is parametrized as = 

Ud cos (TTXi/a) + Mt cos {'KXi/2a + Xt) where Xt is a phase 
factor. The term in eq. (9) is the 4A:F-Umplapp term hav- 
ing the form identical to eq. (6) except the presence of Ud- 
Therefore provided Ud 7^ it can actually produce the 
dimer Mott insulator. On the other hand, in eq. (10), 
gt <x t couples the two phase variables for charge {9p) 
and spin (Oa). In the weakly correlated regime this term 
shows the instability of the system toward the 2fcF CDW 
state. ''^ 

In the strongly correlated regime, on the other hand, 
the competition between the CO insulator {9p = or 
7r/2) and the dimer Mott insulator {9p = 7r/4) arises due 
to the Sfcp-Umklapp term in Hp and the 4A:F-Unrklapp 
term eq. (9), respectively. This is similar to the case in 
§ 2.3 for HiB in eq. (2) with 6d ^ 0, although here, 
Ud, and consequently the presence of the 4fcF-Umklapp 
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Fig. 9. Ground state phase diagram of the one-dimensional ex- 
tended Peierls-Hubbard model at quarter-filling calculated by 
the density matrix renormalization group method, on the plane 
of U/t and V/t for fixed values of 1/K = (dotted line, equiv- 
alent to Fig. 4) and l/K = 1 (filled Unc). CO-SP, DM-SP, and 
CDW phases denote the coexistent state of charge order and 
spin-Peierls lattice tetramization, that of dimer Mott insulator 
and spin-Peierls lattice tetramization, and the weak coupling 
Peierls-Frolich state. [After ref. 91.] 

term itself, should be determined self-consistently. The 
spin degree of freedom in these two states is derived 
to be described by the same phase Hamiltonian Ti,'^ = 
'Ha- — gtUt J dx cosher + '^ci, which is identical to that 
of the SP model. Therefore the two insulating states 
above are both unstable toward the non-magnetic SP 
state with lattice tetramization, ut ^ 0, resulting in co- 
existent states. 

These symmetry broken states are in fact reproduced 
in the numerical calculations. Let us show a ground state 
phase diagram by Kuwabara et al.^^ determined by the 
numerical DMRG method in Fig. 9. The spatial mod- 
ulation of the lattice distortion at each site is deter- 
mined self-consistently and system sizes up to L = 36 
are needed to diminish the finite-size effect. One can see 
that all the states predicted in the phase Hamiltonian 
appear in different regions of parameters {U/t, V/t). For 
the actual pattern of charge density and lattice distortion 
in each state, see ref. 91. 

There is another type of coupling to the lattice de- 
gree of freedom from the Pcierls coupling above: the Hol- 
stein coupling, which describes modulation of the one- 
particle on-site energy at the cost of elastic energy. The 
origin of this coupling can be the electron-molecular vi- 
bration (e-mv) coupling frequently discussed in molec- 
ular systems. This coupling is usually not large such 
that the nature of phenomena due to strong correla- 
tion discussed in this review is qualitatively modified. 
However, recently, exceptionally large deformation of 
molecules themselves coupled to charge disproportiona- 
tion in (EDO-TTF)2PF6 is observed,^'' suggesting large 
e-mv coupling in this material. Another source for the 
Holstein-type coupling recently considered affecting the 
CO state is the influece of the X~ unit implied by experi- 
ments in TM2X.^^''^^ The ferroelectricity observed below 
the CO temperature is interpreted as the combination of 
CO state in the dimerized stack of TM molecules (see 
Fig. 6) and the q = motion of the anions, producing 
large electrical polarization. 




Fig. 10. (Color online) (a) Zigzag ladder and (b) one-dimensional 
chain with next nearest neighbor interactions. These two are 
equivalent with each other as seen in the figure. 

Numerical ED calculations for small cluster sizes up 
to about L — 16 sites including this Holstein coupling 
as well as the Peierls coupling suggests that the former 
type works toward stabilization of the CO state in gen- 
eral. ^^'^^ We note that since the lattice is 3D in nature, 
these couplings would result in finite temperature phase 
transitions,''^ which are not fully investigated yet. Fur- 
thermore, the inclusion of intrinsic dimerization Sd in the 
model might modulate the results, such as a further sta- 
bilization of the dimer Mott insulating state. ^^'^^ Still, 
systematic works which would provide unified pictures of 
the quasi-lD electron-lattice coupled quarter-filled sys- 
tems, together with the interchain coupling discussed in 
§ 2.4, are to be done.^^ 

2. 6 Geometrical frustration 

As introduced in § 2.4, the interchain Coulomb inter- 
action V± added to the ID EHM stabilizes the 2D or 
the 3D CO state with an anti-phase pattern between 
chains. However, in the actual molecular materials the 
manner of coupling between chains is frequently more 
anisotropic. For example, in TM2X, as shown in Fig. 1, 
three kinds of interchain bonds are formed in a zigzag 
way. Then the interchain Coulomb repulsion Vh favors 
the CO pattern along the 6-direction to be anti-phase be- 
tween chains, while Vpi and Vp2 favor the in-phase pat- 
tern. (The observed ferroelectricity in ref. 37 points to 
the latter.) Since these Coulomb energies are expected 
to have similar values, ^^'^^ the two CO patterns com- 
pete with each other. This is an example of geometrical 
frustration. 

Effects of the geometrical frustration are intensively 
studied in localized spin systems^^ defined on lattice 
structures with odd number of sites when rounding along 
the bonds of a cell, e.g., triangular forms as in the above 
example of TM2X. It is known to drive the system to- 
ward destabilization of the classical antiferromagnetic 
state and sometimes even result in a total destruction 
of it, namely, in a spin disordered state. The analogy 
between the classical picture of CO and the spin sys- 
tems noted in § 2.2 naturally leads us to the problem 
of geometrical frustration in CO systems, where the CO 
may be destabilized as well. This was first considered by 
Anderson^^ in the context of the "classical" 3D CO sys- 
tem, magnetite Fe304, described as an Ising model on 
the phyloclore lattice which is a typical frustrated lattice 
structure. 

Treating the zigzag coupling in the 2D model for 
TM2Ar with full anisotropy is a problem yet to be 
solved. A similar situation will be discussed in § 3.3, re- 
lated to the CO in quasi-2D 9-KI2X. Instead, here we 
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consider a more simple system but having many com- 
mon aspects: a double-chain system with zigzag cou- 
plings between them, i.e., a zigzag ladder, shown in 
Fig. 10. As seen in Fig. 10, it is equivalent to a ID 
model with next-nearest-neighbor interactions. Then the 
EHM on this chain structure can be written as Hib 
in eq. (2) with (5d = added by next-nearest terms, 
*2 J2iai4a^i+2cr + h-C.) -I- V2 J^i nini+2- In the following 
we write the nearest neighbor terms as ti and Vi. 

This model at quarter-filling has been studied as an 
effective model for two different material systems. One 
is for TMTSF2X, with t2 = 0, extending the range of 
Coulomb interaction terms from the previous works. ^"^'^^ 
Another is for the Cu-0 subunit of a transition metal ox- 
ide PrBa2Cu408, where the effective model for the Cu 
sites has a zigzag ladder structure with ti <ti ^2, and 
Vi ~ V2V2 suggested by the Cu-Cu distances. 1°"' 
For ti = or ^2 = 0, the model at J7 = 00 is exactly 
mapped onto XXZ models when one follows the mapping 
of Ovchinnikov explained in § 2.2, therefore the analogy 
between the frustrated localized spin system and the CO 
system is straightforward here again. 

For the t2 = case, Emery and Noguera^^^ actually 
studied the ID XXZ model with next nearest Jz term, 
while Yoshioka et alJ^ studied the finite- C/ case by the 
method described in § 2.3. In these studies the interaction 
terms are treated perpurbatively from the weak-coupling 
regime. They have discussed the competition between 
two different CO patterns favored respectively by Vi and 
V2. The one favored by V2 having wave vector 2A:f (^f 
defined for t2 — 0) shows 2k-p SDW correlation devel- 
opped. The co-existent state of 2fcF SDW and 2/cf CDW 
was first proposed based on the MF approximation,^^ 
to discuss such state observed in TMTSF2Ar.^° We note 
that this "2A:f" CDW is different from the Peierls-Frolich 
state due to electron-lattice coupling, but it is originated 
from the strong correlation, V2, i.e., it is indeed a CO 
state. A possible insulating bond ordered wave phase is 
also found in the above weak-coupling approaches, in the 
competing region. 

Numerical studies on such models on the zigzag ladder 
have found a different aspect at large X^j's where the 
weak coupling approaches above may break down: a wide 
region of the metallic TLL phase between the two CO 
phases. A phase diagram by the numerical DMRG 
method is shown in Fig. 11,^°^ where TLL phase is widely 
realized in the Vi ~ 2V2 region, with Kp close to 1/4 
suggesting the closeness to CO. This is a consequence of 
the geometrical frustration, and in such metallic phase 
with both Vi and V2 having large values, the existence 
of large fluctuations of both CO states is expected. 

Recently, it is pointed out that the CO state in DL 
DCNQl2Ag is affected by similar geometrical frustra- 
tion. This is due to a peculiar situation where the inter- 
chain configuration along each plaquette in the a-h plane 
(see Fig. 1) has a spiral symmetry which is not com- 
patible with the twofold periodicity generated by CO.^"^ 
Such frustration due to the spiral symmetry is also found 
to exist in the process of the metal-insulator transition 
in DMe-DCNQl2Cu,^°^ where the 7r-electron couples to 
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Fig. 11. Ground state phase diagram of the extended Hubbard 
model at quarter-filling on the zigzag ladder, determined by the 
density matrix renormalization method. ^'''^ is fixed at 10 

while (2 is set to 0. Sandwitched by the two insulating CO 
phases with different charge patterns, the Tomonaga-Luttinger 
liquid (TLL) metallic phase is extended toward the frustrated 
line Vi = 2V2. The contour lines indicate the TLL parameter 
Kp. [Reprinted figure with permission from S. Ejima et al, Phys. 
Rev. B 72 (2005) 033101. Copyright (2005) by the American 
Physical Society.] 

the Cu c?-electron producing a three-fold periodic struc- 
ture^'^'' again not compatible with the spiral symmetry, 
which is beyond the scope of this review. 

3. Quasi-Two-Dimensional Systems 

Many theoretical works on CO states in 2D models 
near half-filling have been performed, such as in the 
Hubbard model and in the t-J model on the square lat- 
tice, mainly devoted to the high-Tc SC cuprates and re- 
lated transition metal oxides. "'^"^ However, studies on the 
quarter-filled case are mostly triggered by the recent ex- 
perimental findings of CO in quasi-2D molecular conduc- 
tors. To the authors' knowledge, only a few have been 
done beforehand. Ohta et al}°^ studied the 2D EHM 
on the square lattice for different fillings including one 
quarter, in fact motivated by the cuprates, while a more 
specific case was studied^ to explain the CO transition 
in NaV205 (see § 4). 

Since CO has been found in many quasi-2D A2B sys- 
tems such as in ET2X, theoretical works have been per- 
formed on the 2D EHM adopting the full anisotropy (or 
a part of it) of the materials. On the other hand, the 
quarter-filled EHM on the square lattice has been studied 
as well from a rather general viewpoint. Such 2D models 
are difficult to analyze theoretically in a controlled way, 
and intensive efforts are now in progress. We will review 
the present status in this section. 

3.1 Charge ordering in quasi-two-dimensional systems 

Direct observation of the CO transition in quasi-2D 
A2B compounds was first in 6l-ET2RbZn(SCN)4,''' 
and next in a-ET2l3,"^ both found in ^^C-NMR exper- 
iments. Now many quasi-2D A2B materials, including 
non-ET based compounds, are recognized to show CO.^^ 
Their crystal structures are composed of alternating lay- 
ers of A molecules and B units, and quasi-2D electronic 
states are realized; tij along the interlayer direction are 
small due to the insulating B layers. These materials 
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Fig. 12. (Color online) Schematic representation of the 2D layers 
in different polytypes of quasi-2D A2B materials. The indices 
are not only for tij but also for different values of Vij, although 
the degree of anisotropy can be different (see text). In the re- 
type and the A-type structures the thick A bonds are those with 
the transfer integral Ia considerably larger than the others, i.e., 
dimers. [By courtesy of C. Hotta.] 

have a variety in their 2D arrangements of A molecules, 
sometimes even for the same chemical formula, i.e., poly- 
types, classified by greek characters, a, /?, k, 0, A, and so 
on.^^ Some representative examples are shown in Fig. 12, 
for which the actual values of transfer integrals calculated 
by the extended Hiickel method we refer to, e.g. refs. 8 
and 45. Because of such anisotropy in the lattice struc- 
tures, the EHM we should treat, TYehm in eq. (1), would 
be complicated. As was mentioned in § 1, we emphasize 
that such diversity of lattice structures here is one of the 
important characteristics producing the rich variety of 
properties in molecular conductors. CO is in fact very 
much affected by such geometry of the lattice as we will 
see in the following subsections. 

A theoretical approach to handle this diversity was 
pursued by Kino and Fukuyama,^^ who considered Hub- 
bard models {Vij = 0) for ET2X, taking into account the 
full anisotropy in Uj for each polytype, and applied MP 
approximation to the on-site Coulomb interaction term 
U. By this they have provided a way to search for rela- 
tions between the crystal structure and electronic prop- 
erties, and also for generic understandings even among 
different polytypes. Since the polytypes show a va- 
riety, different ways of such systematic views have been 
discussed by different authors since then.^''*''^' ^^"^ One way 
is to consider an anisotropic triangular lattice, re- 
fined by Hotta^^^ as follows. 

In Fig. 12, one can see that 9, a, and /3-type struc- 
tures have the triangular lattice as basis, "colored" 
with anisotropy. Since the charge transfer is realized as 
[ET2]^A'^, the one-particle HOMO bands as a whole are 
quarter- filled in terms of holes. Effective models for the 



materials with these symmetries would be the quarter- 
filled EHM on an anisotropic triangular lattice, with Vij 
having the same anisotropy as Uj. For example, for the 
9-type structure with highest symmetry among the poly- 
types, 

ne=~tp (cLcj<T + h.c.) -tc (cLcj<, + h.c.) 

+ UY^iT^il + X! ^^"^3 + ^ riiTlj, (11) 

where {ij)p and {ij)c denotes the site pairs i and j along 
the p bonds and the c bonds, respectively. In contrast, 
K and A-type stuctures consist of dimers, i.e., pairs of 
molecules connected by Uj considerably larger than the 
others, indicated as the A bonds in Fig. 12; if one consid- 
ers the dimers as units, these structures are topologically 
equivalent to 9 and /3-type structures, respectively. Half- 
filled models on the anisotropic triangular lattice can be 
effective models for such systems with large dimeriza- 
tion. Then the on-dimer Coulomb repulsion J/dimer can 
lead the system toward the dimer Mott insulator (see 
refs. 8, 17-20). 

The CO states are stabilized in the former quarter- 
filled models. In Kino and Fukuyama's work on the Hub- 
bard model for Q;-ET2l3,^^^ an insulating state accompa- 
nied with charge disproportionation between stacks (the 
vertical stripe type CO pattern; see § 3.3) is found. How- 
ever, the charge pattern there is not consistent with later 
experiments. Moreover in a MF study on the Hubbard 
model for the other CO material 6i-ET2RbZn(SCN)4, CO 
state is not even stabilized. Following these works, Vij 
added to the Hubbard model retaining the full anisotropy 
in tij, i.e., the EHM for ET2X on the anisotropic triangu- 
lar lattice, has been introduced by Seo,^^* who discussed 
the CO states there based on the MF approximation. As 
in the quasi- ID systems discussed in § 2, the intersite 
Coulomb interaction is also crucial for CO in the quasi- 
2D materials. 

The values of Vij are in fact estimated to be appre- 
ciable compared to U. As in the quasi-lD compounds in 
§ 2, U is believed to be of the order of 1 eV in the FT 
compounds as well, while, again, the actual values for 
Vij/U estimated as 0.2 ^ 0.7^^^"^^^ are too ambiguous 
for quantitative arguments. Typical values for \tij \ in the 
FT materials estimated by the extended Hiickel method 
range around 0.1 ~ 0.25 eV.*'"*^ Therefore these ma- 
terials are indeed strongly correlated systems. We note 
again that the degree of anisotropy in tij and that in Vij 
does not correspond in a one-to-one manner. For example 
in many members of 6'-ET2A' described by eq. (11), \tc\ 
is much smaller than |ip|,^'* despite of their close inter- 
molecular distance resulting in ~ Vp}^^ 

The anisotropic triangular lattice is somewhat com- 
plicated and therefore it is not only tough to handle it 
theoretically but also sometimes difficult to extract ex- 
plicit key parameters for the physics therein. One natu- 
ral way of thinking, as was pointed out by McKenzie et 
al.,^^^ is to study the quarter-filled EHM on the square 
lattice, which is much simpler but would still show essen- 
tial features of the CO transition in 2D, and even of the 
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expected SC state^^'^ (§ 3.4). In this model, there is one 
kind of transfer integral t and also inter-site Coulomb 
repulsion V, along nearest-neighbor sites both in x and 
y directions. Such studies on the square lattice can shed 
light on the properties of the anisotropic triangular lat- 
tice EHM, as the former is one of the limiting cases of 
the latter. For example, in the 0-type structure if we set 
tc = 0, Vc = and t = tp, V = Vp, Tig in eq. (11) be- 
comes equivalent to the square lattice EHM. Neglecting 
tc is in fact a good approximation for whereas 
Vc cannot be neglected, as mentioned above, while this 
square lattice case might be applied to some /3"-type 
compounds. Note that electron-hole symmetry holds 
in the square lattice EHM while tc breaks it, so relative 
signs between tij become distinct. We will explain stud- 
ies on the square lattice EHM first (§ 3.2) and then the 
anisotropic triangular lattice case next (§ 3.3), although 
the research has been developped more or less in a par- 
allel way. 

3.2 Extended Hubbard model on square lattice 

The existence of a checkerboard type CO with wave 
vector q={-K / a,T: / a) on the quarter-filled square lattice 
EHM, as shown in Fig. 13, has been first demonstrated 
by Ohta et aL,^°^ for U = 8t and V = 3t. It is based on 
calculations of the equal-time charge correlation function 
C{q) = '^ij{ninj)e'^'i'^^^ , where L is the total num- 
ber of sites in the square shaped clusters and Rij is the 
vector connecting site i and j, by use of the Lanczos ED 
technique with cluster sizes up to L = 16. This pattern 
can be understood naturally as the extention of the ID 
case (Fig. 2(a)) from the idea of "Wigner crystal on lat- 
tice" . 

They have also derived an effective spin model in such 
CO state by the fourth order perturbative expansion 
from strong coupling, t <C {U,V). It is an antiferromag- 
netic Heisenberg model on the square lattice, rotated 45 
degrees from the original lattice (see Fig. 13), with a 
nearest neighbor exchange coupling, 



J 
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Note that even at U /t — 00, novel ring exchange pro- 
cesses break the zeroth order spin degeneracy resulting in 
a finite J = At'^ /9V^ Since the next-nearest-neighbor 
coupling J' is estimated to be much smaller, an antifer- 
romagnetic spin order as shown in Fig. 13 is expected 
at the ground state since the square lattice Heisenberg 
model shows Neel order. 

A systematic numerical Lanczos ED study was demon- 
strated by Calandra et aL,^^^ for cluster sizes up to 
L = 20. In this case, a powerful method to determine 
from small systems whether the bulk system is metallic 
or insulating is to evaluate the Drude weight £). 126-128 
It is given by 

(0|i^|0) |(n|j,|0)|2 
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27re2 
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L ^ En — Eq 
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Fig. 13. (Color online) Checkerboard charge ordered insulating 
ground state of the square lattice extended Hubbard model at 
quarter filling in large U/t and V/t. An antiferromagnetic inter- 
action J occurs between spins along the diagonals shown in the 
figure. 




where Eq and £"„ denote the ground state (|0)) and ex- 
cited state (|n)) energies of the system, respectively.^^ K 



Fig. 14. The Drude weight, D, for the quarter-filled square lattice 
extended Hubbard model as a function of V/t, for L = 16 and 
various values of U/t (continous and dotted lines) , and for L = 
20 and U = lOt (dashed line).^'^^ The arrow in the horizontal 
axis marks the onset of checkerboard charge ordering for U = 
lot. The inset shows the finite-size scaling of D as a function 
of 1/Lx for different values of V/t with U = lOt. The metal- 
insulator transition occurs at V^^ ~ 2.2t. [Reprinted figure with 
permission from M. Calandra et ah, Phys. Rev. B 66 (2002) 
195102. Copyright (2002) by the American Physical Society.] 

is the kinetic energy operator and jx is the current oper- 
ator in the x direction at zero wavevector {q — 0). The 
occurrence of an insulating phase is marked by the expo- 
nential vanishing of D with the linear size of the system 
Lx = VL.^^^'^^^ In Fig. 14 the results are shown, where 
D is plotted as a function of V/ 1 for different U/t. As 
V/t is increased, D decreases until it eventually vanishes. 
For U/t = 10 the critical value for the metal-insulator 
transition is estimated as V^^/t ~ 2.2. We should note 
that this might not be conclusive as an accurate value 
as the finite size scaling shown in the inset of Fig. 14 
is only for the three clusters L=8, 16, and 20, due to 
computorwise limit. Nevertheless, for V > V^^, D dis- 
plays an exponential dependence with 1 / as expected 
for an insulator, while for V < V^^, in contrast, 
it weakly depends on 1 / extrapolated to a finite value 
in the thermodynamic limit, consistent with a metallic 
state. On the other hand, the occurrence of CO has been 
investigated by computing C{q) systematically. The ex- 
trapolation suggests V^'^/t ~ 1.6, where C(7r/a, 7r/a) 
starts to show a peak as V is increased, which possibly 
overestimate its value due to finite size effects. When U 
decreases the critical value V^'~' increases, similarly as in 
the ID case (see Fig. 4). 
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Since the estimates above provide > V^^'--' , a pos- 
sible CO metallic phase, not found in the ID EHM, is 
suggested. In fact, a recent variational Monte-Carlo cal- 
culation-'^^^ (see next subsection) have confirmed its ex- 
istence, with a CO metallic phase realized between the 
uniform metallic and the CO insulating phases, i.e., at 
V^r° < V" < V^^. We note that another work using 
projector correlated method^^° shows that such feature 
is also seen at finite temperatures. This is in clear con- 
trast with the purely ID EHM discussed in § 2, where 
the quantum fluctuation destroys the finite-T CO phase 
transition. More theoretical efforts are needed, to explore 
the overall phase diagram for this model, in the (U, V, T) 
parameter space, for the critical values of CO, metal- 
insulator, and antiferromagnetic phase transitions. 

Although the analysis above is restricted to numerics, 
the U/t = 00 limit has been investigated analytically by 
McKenzie et al}"^^ In this limit the square lattice EHM 
reads after projecting out the doubly occupied sites as 

Htv = -tJ2'^ + h.c^j P + VY, rurij, (14) 

where V projects out the doubly occupied configurations 
due to their prohibition at U/t = oo. We refer to this 
model as the t-V model; this is the generalization of the 
U/t — oo limit of the ID case discussed in § 2.2, whereas 
the charge and spin degree of freedoms cannot be decou- 
pled in 2D. 

They have introduced an SU(iV) generalization of this 
model in which the spin index runs from 1 to N, and 
carried out a slave boson theory, which have been ap- 
plied to various models for strongly correlated electron 
systems by different authors. The electron creation 
operator is replaced by = fja^it where the spinless 
charged boson operator bi is introduced to keep track of 
the empty sites, and //^ is a fermion operator carrying 
spin. In order to preserve the anticommutation relation 
for the electrons the new operators must satisfy the local 
constraint fl^fia + b\bi = N/2. Hereafter, whenever a 
repeated a index appears a sum from 1 to A'' is assumed. 
Note that the original EHM has A'^ = 2. 

In the coherent state path integral formulation^^^ the 
Lagrangian at imaginary time r is given by 

i 

{ij) (u> 

+ T.iWlfia + blbi-N/2), (15) 

i 

where /x is the chemical potential and A, is a static La- 
grange multiplier enforcing the local constraint above. 
We have used the fact that cl^Cia = fl^fia- 

Following Kotiar and Liu's work on the square lattice 
Hubbard model, ^"'^ it is convenient to choose the radial 
gauge to avoid possible infrared divergences where the 
boson amplitude becomes a real number, rj = |6j|, and 
Aj becomes a dynamical bosonic field: Xi{T). By using the 
relation fj^fia = N/2—h\bi to replace one pair of fermion 
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operators in the V term and integrating out the other 
fermionic degree of freedoms, we are left with an effective 
Lagrangian which is quadratic in the boson fields. 

The MP solution of the effective bosonic model is 
obtained by assuming that the boson fields are spa- 
tially homogeneous and time independent: rj(r) = b and 
iXiir) = A. This treatment is exact in the A' — > oo 
limit, therefore is a base for the expansion in powers of 
1/A'.^^^'^^^ The resulting MF free energy is 

where /3 = l/{kBT) and ujn is the fermion Matsub- 
ara frequency. The MF eigenenergy is given by — 
-2tb'^ (cos kx+cosky)/N + X- + AVn/N. Minimization 
of F^^ with respect to b and A leads to 6^ = N/2-n, A = 
2t^f^ /(efc)(cosfcj. + cos ky + 4V). and /i is adjusted to 
give the correct electron filling. This form, eq. (16), indi- 
cates that the MF solution for the large- A' generalized t- 
V model describes just a rcnormalizcd Fermi liquid, anal- 
ogous to the case of the Hubbard model. ^'^^ In the case of 
quarter-filling and N = 2 (which is somewhat artificial 
as one implicitly consider A to be large) , the bandwidth 
is reduced to half its bare value as 6^ = 1/2 and the band 
is shifted from its bare position by A. Namely, the mass is 
enhanced as twice the bare value, i.e., m*/m = 1/5^ = 2. 
The overall effect of the nearest-neighbour Coulomb in- 
teraction, V, reduces to a constant shift in the chemical 
potential. 

The leading 1/A^ corrections modify the MF solution 
so that when V/t is increased the Fermi liquid phase re- 
sults in an instability toward the checkerboard type CO 
state. This is seen when writing the boson fields in terms 
of the static MF solution (b, A) and the dynamic fluctu- 
ating parts: rj(r) = b + bSri^r), and iXii^r) = X + iSXi^r). 
The resulting effective action, to the second order in the 
boson fields, is 5 = + 5^^^ , where the second term 
due to the fluctuations in the boson fields is written in 
the form as 




where Vn is the boson Matsubara frequency; for the 
explicit form of the elements in r(q, i/n) see ref. 122. 
We note that F^a is the Lindhard function describ- 
ing density-density fluctuations in the renormalized 
band. The propagators of the boson fields D{q, Vn) = 
T~^{q, Vn) are of the order 0{1/N), as they should. The 
condition for an instability of the Fermi liquid phase is 
when the quadratic form (17) becomes negative at some 
q for 1/ = 0, so that fluctuations in the charge den- 
sity will decrease the free energy. This is determined by 
detr < at q = {n/a, tt/o) for quarter filling, which pro- 
vides {V/t)cT = 0.69 for the critical value from the Fermi 
liquid to the CO state for this largc-A t-V model. We 
note that the inclusion of diagonal hopping t' (equivalent 
to tc in eq. (11)) is easily incorporated in this scheme, 
which does not modify the results qualitatively. 
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Fig. 15. (Color online) Charge ordered states in the quarter-filled 
extended Hubbard models on anisotropic triangular lattices: (a) 
vertical, (b) diagonal, and (c) horizontal stripe type states, and 
(d) the threefold state. Sites with/ without colored circle rep- 
resent the charge rich/poor sites. The solid and dotted bonds 
represent the "p" (pi ~ p4) and "c" (cl ~ c3) bonds in the 9(or 
^'d)-typo and the o-type structures in Fig. 12. 

3.3 Extended Hubbard model on anisotropic triangular 
lattice 

As we have seen in the previous subsection, the inter- 
site Coulomb repulsion can drive the system toward CO 
in 2D systems as well. In the strongly correlated regime, 
the spins are described by an effective 2D Heisenberg 
model with spins located on the charge rich sites. Such 
aspects are common with the ID cases discussed in § 2, 
where CO results in magnetic properties of ID Heisen- 
berg chains. These knowledges are very useful in treating 
the more anisotropic cases introduced in § 3.1, with full 
anisotropy taken into account. There, calculations be- 
yond MF are still in progress, but if we deduce from the 
MF results together with the picture above, qualitative 
understandings of the experiments can be obtained^^ as 
follows. 

MF calculations on the anisotropic triangular lattice 
EHM were first performed by Seo^^^ for different mate- 
rials with 9 and a-type structures. The U and the Vij 
terms are treated within the standard MF treatment as 
) -f {nia)n-ja' - {nia){nja'), by allow- 
ing large unit cell sizes to consider various CO states 
with several possible spin orders for each of them. In 
fact many self-consistent MF solutions are obtained and 
their ground state energies are compared. It is found that 
CO states with lowest energy show a variety as shown in 
Fig. 15 depending on the parameters, which is in con- 
trast to the square lattice case where the checkerboard 
pattern is stable. Mainly the "stripe" type CO states are 
discussed in ref. 118, shown in Fig. 15(a) ~ (c), where 
the carrier densities on the charge rich and the charge 
poor sites are (in some cases approximately) 1/2 + d and 
1/2 — S, respectively. These can be considered as exten- 
sions of the Wigner crystal- type CO, which can lead to 
insulators in the strongly correlated regime. In general, 
the vertical stripe type pattern is favored by the Coulomb 
repulsion along the solid bonds in Fig. 15, while that 
along the dotted bonds stabilizes the diagonal and hori- 
zontal patterns. 

Let us start with the MF results for the 6'-type struc- 



(a) 6lj-ET2RbZn(SCN)4 (b) a-ET2l3 

mm mH 

Fig. 16. (Color online) Schematic drawing of mean-field solu- 
tions^^** for the horizontal stripe-type charge ordered states in 
extended Hubbard model on anisotropic triangular lattices for 
(a) ed-ET2RbZn(SCN)4 and (b) 0-ET3I3. Different "p" bonds 
are represented by different lines following Fig. 12. 

ture, eq. (11), which has the highest symmetry among 
the polytypes. As noted in § 3.1, nearly isotropic Vc — Vp 
is expected, which gives rise to close competition be- 
tween the different CO patterns. MF calculations^^^ for 
eq. (11) at T = considering the stripe type CO patterns 
as in Fig. 15(a) ~ (c), have been performed for tp = —0.1 
eV, tc = 0.01 eV, taken from typical extended Hiickel 
parameters for 6'-ET2MM'(SCN)4 (M: Rb, Cs, Tl, etc., 
and M': Zn, Co, etc.),^"* and U = 0.7 eV, while varying 
Vp and Vc near 0.3 eV. The results show that either ver- 
tical or diagonal stripe type CO solution is stabilized for 
^c'^Vp and Vc>Vp, respectively, whereas the horizon- 
tal pattern has rather higher MF energies. In the actual 
0-type compounds which sustain this structure down to 
low temperatures, either vertical^'^^ or diagonal^^"* pat- 
tern of CO is indicated by optical measurements, consis- 
tent with the MF results. We note that the vertical stripe 
CO state here is equivalent to the checkerboard pattern 
in the = and = square lattice EHM discussed 
in § 3.2. 

In the region of Vc — Vp, first in a study setting 
tij = O^'^^ and recently in a MF study, ^'^^ a solution 
with a larger periodicity called as the "threefold" state 
(Fig. 15 (d)), which was not considered in ref. 118, is 
found to have lower energy than the stripe-type CO 
states. In this state, the carrier densities on the charge 
rich and the charge poor sites are 1/2 + 25 and 1/2 — (5, 
respectively (the tij = study^^^ provides the extreme 
case of (5 = 1/4, namely, carrier densities of 1 and 1/4). In 
this state, the carriers can avoid the occupancy in both 
neighboring pairs, Vc and Vp, therefore it is stabilized 
ai Vc — Vp, compared with the stripe- type CO states 
which cost loss in either one of these Coulomb repul- 
sion terms. Such a long period state in fact appears as a 
short range order as found in X-ray scattering measure- 
ments^^^' ^'^^ as diffusive rods (no interlayer coherence), 
in the conductive states of several members of 0-ET2X. 
However these are observed at different 2D wave vectors 
than in the studies above, which may be resolved by in- 
troducing longer range of Coulomb interactions in the 
EHM.-'^'^^ In any case this threefold state is metallic since 
it is incompatible with the periodicity for quarter-filling, 
in contrast with the stripe-type CO states which is insu- 
lating in the strongly correlated regime, therefore these 
instabilites should be distinguished. 

Many ^-type compounds exhibiting CO undergo a 
first order phase transition accompanying a structural 
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transformation to a slightly dimcrizcd structure, typi- 
cally seen in several members of 6'-ET2MM'(SCN)4:2* 
the ^d-type structure, with a modulated molecular net- 
work as shown in Fig. 12. MF calculations^"'^®' "'^^^ under 
such structure with full anisotropy of tij taken into ac- 
count, but approximating Vij to be two kinds, Vp and 
Vc, as in the calculation for the 0-typc structure, have 
been performed. The results suggest that a horizontal 
stripe- type solution as shown in Fig. 16(a) can have the 
lowest energy for 14 > V^. This is the pattern observed 
in these compounds, such as in ^-ET2RbZn(SCN)4 below 
the metal-insluator transition temperature Tmi-^^® Such 
calculations have been performed by directly adopting 
the low temperature structure for simplicity, while the 
importance of the coupling to the lattice have been im- 
plied.^' It is interesting to start from the high tem- 
perature ^-type structure and see whether the structural 
phase transition to ^d-typc together with the CO transi- 
tion could be reproduced, which awaits future studies. 

As for the a-type structure adopting the values of tij 
for a-ET2l3, MF calculations^^* show that a horizontal 
stripe solution (Fig. 16(b)) is more stable than the other 
solutions in the Vc>Vp region, in contrast to the verti- 
cal stripe type solution found by Kino and Fukuyama^^^ 
stabilized in the small Vij region. This horizontal stripe 
pattern is in fact in accordance with the experiments be- 
low the metal-insulator CO phase transition temperature 
in this material. It is noteworthy that, in both this salt 
and in the 9d-type ET2X above, the MF calculations and 
the experiments agree with each other to a degree of the 
actual molecules having rich and poor carrier density in 
their rather complicated unit cells. 

Recently, Kobayashi et al}'^^'^^^ have extended such 
MF calculations on this compound, a-ET2l3, to param- 
eters under pressure, where a so-called narrow-gap semi- 
conducting state is found by experiments.^^*' The re- 
sults show successive transitions from the CO insulat- 
ing to the CO metallic phase, and then to a phase with 
a peculiar "zero-gap" state, as a function of pressure. 
This zero-gap state is characterized by a dispersion of 
anisotropic Dirac fermions at a wave vector fe°, which 
crosses at the Fermi energy. It is noteworthy that fc° is an 
incommensurate wave vc^etor due to an "accidentia! de- 
generacy" naturally emerging in the band structure of 
this compound, but not related with the symmetry of the 
lattice striicture. There are attempts^^"* to explain the 
anomalous transport properties observed at high pres- 
sures^**' based on such framework. 

The magnetic properties under the stripe-type CO 
states can be deduced from such MF results together 
with the picture mentioned above that we have learned 
from the ID cases and the 2D square lattice case. The 
spin degree of freedom is expected to be described 
by (quasi-) ID Heisenberg systems with spins on the 
charge rich sites along the stripes, then quantum fluc- 
tuation may destroy the antiferromagnetic spin order 
found in MF solutions (see Fig. 16). In fact, such be- 
havior is typically observed in the measurements on 6- 
ET2RbZn(SCN)4 and a-ET2l3 at ambient pressure, both 
showing horizontal stripe type CO but absolutely differ- 
ent magnetic behavior. In the former, the CO is real- 



ized in the ^^-type structure, where the MF results show 
that the bonds between charge rich sites are all equiv- 
alent along one kind of bond, p4 in Fig. 12 (the double 
solid bonds in Fig. 16(a)), resulting in a uniform Heisen- 
berg coupling Jp4 between spins along these stripes. On 
the other hand, in the latter compound, the horizontal 
stripe pattern results in an alternation of bonds along 
the charge rich sites, p2 and p3 in Fig. 12 (the thick 
solid and long dashed bonds in Fig. 16(b)), therefore 
spin singlet formation due to alternating Jp2 and Jps 
can be expected. These explain the difference seen in 
the magnetic susceptibility data: a Bonner-Fischer like 
low- dimensional localized spin behavior in the former^* 
and a prominent spin-gap behavior in the latter.^*'' We 
note that in ^-ET2RbZn(SCN)4, at lower temperature 
another phase transition takes place: a spin singlet for- 
mation due to the spin-Peierls mechanism along the 
stripes.^* This is another experimental fact indicating 
the validity of describing the magnetic properties under 
stripe-type CO states as ID Heisenberg models. 

Although such procedure above is helpful in qualita- 
tively understand the experiments, MF calculations can- 
not accurately describe situations where quantum fluc- 
tuations play crucial roles as mentioned in § 2. This is 
indeed the case for the 0-type structure near Vc — Vp, 
as this is a situation under geometrical frustration intro- 
duced in § 2.6. In fact, computations including quantum 
fluctuation by Merino et al}'^^ showed that the competi- 
tion between different CO patterns can lead to a "quan- 
tum melting" of CO, which is analogous to the case of 
the zigzag ladder model mentioned in § 2.6. The ground 
state phase diagram on the Vc-Vp plane for fixed U = lOt 
is shown in Fig. 17, which is based on the Drude weight 
calculations on a L = 16 site cluster within the Lanc- 
zos ED technique explained in the previous subsection. 
There, a large metallic phase is stabilized in the region 
of Vc ~ as in the ID case we have seen in Fig. 11. 
The two different stripe-type CO states are stabilized in 
the less frustrated regions Vc ^ Vp and Vc <^Vp, which 
are, in contrast, consistent with the MF studies explained 
above. 

Recently another study for this EHM on the ^^-type 
structure using the variational Monte Carlo method 
has been performed for much larger systems. In this 
method, variational states with CO patterns together 
with correlation effects are assumed and the variational 
energy is calculated numerically by Monte Carlo sam- 
pling. The features in the ED calculations above are 
reproduced, such as the robust metallic phase along the 
Vp ^ Vc line. In this metallic phase, the threefold pe- 
riodic modulation of charge density is found as in the 
MF study^'"' noted above, which, however, did not fit 
the L = 16 cluster used in the ED study. In the large 
Vp '^Vc region, the threefold modulation is rather strong, 
which becomes weaker for small Vp "^Vc- There, different 
metallic states with stripe- type modulation, i.e., stripe- 
type CO metallic states, have very close variational en- 
ergies. This indicates that there is large charge fluctua- 
^jQj^i04, 148 inciu(iing the threefold state due to the geo- 
metrical frustration. 

Combining the above results of different theoretical 
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Fig. 17. Lanczos exact diagonalization ground state phase dia- 
gram of the quarter-filled extended Hubbard model on the 6-type 
lattice structure for {/ = lOt.^''* The figure is based on "square 
lattice" representation: V and V is respectively equavalent to 
and Vc in eq. (11), and tp = t, tc = 0. The boundaries are 
extracted from the Drude weight calculation on a L = 16 site 
cluster. The label M stands for metal, and COl for charge or- 
dered insulator. [Reprinted figure with permission from J. Merino 
et al: Phys. Rev. B 71 (2005) 125111. Copyright (2005) by the 
American Physical Society.] 

methods, one can consider that the horizontal stripe 
CO states found in the actual 6-F,T2X salts arc real- 
ized by relaxing the frustration effect through the struc- 
tural phase transition. In clear contrast, experiments 
show that members of d-'ET2X which does not show 
such structural phase transition are rather conductive, 
but show a gradual increase of resistivity at low temper- 
atures.^'' There are indications of a glassy CO state in 
^^C-NMR measurements, which seems to be related 
with such large charge fluctuation. The additional disor- 
der effect may be responsible for such behavior, which 
awaits to be understood. 

3.4 Superconductivity 

Molecular conductors exhibit a variety of SC states 
as found in quasi-lD Bechgaard salts and quasi-2D ET 
salts.® The issue of their mechanism, particularly the role 
of spin and charge fluctuations has attracted much atten- 
tion. When an SDW or an antiferromagnetic phase 
is located next to the SC phase, it has been often as- 
serted that the spin fluctuation is its origin. As an ex- 
ample, K-ET2X shows SC next to an antiferromagnetic 
Mott insulating phase duo to the dimerization.^" Theo- 
retical calculations for the half-filled Hubbard model on 
the anisotropic triangular lattice for the K-type structure 
show that the SC state here has the dx^-y2 symmcitry me- 
diated by antiferromagnetic spin fluctuation.^^^"^^'^ Sim- 
ilarity to the high-Tc cuprates has been discussed, and 
in fact the rf-wave state is stable when the lattice struc- 
ture is continuously varied from the triangular lattice to 
the square lattice. ^^^'^^^ 

On the other hand, recent experiments on 2D A2B 
salts suggest the existence of a SC phase in the vicinity 
of the CO phase, such as in the unified phase diagram 
of 6'-typc compounds^'' and SC is actually observed in 
6i-(DIETS)2Au(CN)4 under uniaxial pressure.i^'^ In a- 
ET2I3 under uniaxial pressure, SC is even implied to 
show up in the presence of CO.^" These motivated the- 
oretical studies to investigate possibilities of charge fluc- 



tuation as an origin of the SC state, based on the 2D 
EHM on the square lattice as well as on the anisotropic 
triangular lattice, which we will review in this subsection. 

SC mediated by charge fluctuation was discussed flrst 
by Scalapino et al}^^ for the 3D EHM on the cu- 
bic lattice. Relevance of this mechanism to the molec- 
ular conductors was pointed out flrst by Merino and 
McKenzie,^^^ in the quarter-fllled square lattice EHM 
by extending the slave-boson theory for the large- A?' t-V 
model^^^ introduced in § 3.2. Later several authors^^^' 
have extended the weak coupling approach by Scalapino 
et al. (see below) to the square lattice case as well. Their 
results all show that in the 2D case charge fluctuation 
due to the nearest-neighbor Coulomb repulsion V can 
induce a singlet SC state with dxy-wave symmetry. 

In the slave-boson theory for the large-A^ t-V model, 
the effective interaction 14ff(g = k — k'), acting between 
two quasiparticles with momentum k and —k which scat- 
ter to fe' and — fe', can be calculated from the 1/A^ fluctu- 
ations around the MF solution at N 00.^^^'^®^ As the 
ratio V/t is increased, the potential varies its shape devel- 
oping singularities at (±7r/a, ±7r/a) toward the critical 
value {V/t)cr for the occurance of checkerboard CO, at 
zero temperature (see § 3.2). This leads to an attraction 
in the channel with rfj;j,-symmetry, and the SC instabil- 
ity is estimated by the coupling averaged over the Fermi 
surface turning from positive to negative by decreasing 
temperature.^®^ We note that the Fermi surface here has 
poor nesting, then the pairing mechanism is induced by 
the intersite Coulomb repulsion V. 

The resulting T-V phase diagram is shown in Fig. 18, 
where the d^y-wave SC state is next to the CO phase dis- 
cussed in § 3.2, near the critical point {V/t)ci- We note 
that a re-entrant behaviour for the CO phase transition 
temperature is found, which is also seen in the EHM 
in inflnite dimension using the dynamical MF theory. 
From the strong coupling viewpoint, the spin fluctua- 
tion associated with the fourth order exchange process 
in eq. (12) couples the next-nearest neighbor sites anti- 
ferromagnetically (see Fig. 13). This spin interaction, in 
turn, cooperates with the charge fluctuation in stabiliz- 
ing the dxy-wave SC state close to the CO transition. 

Search for a SC state near CO in the anisotropic tri- 
angular lattice EHM has been pursued by Tanaka et 

164 They used the RPA treatment adopting tij for 
6l-(DIETS)2[Au(CN)4] as tjtp = 0.4, and assumed an 
isotropic inter-site Coulomb repulsion Vc = Vp = V, in 
eq. (11). In this approximation, the ladder- type diagrams 
containing V are neglected, but a recent numerical calcu- 
lation in the fluctuation exchange approximation retain- 
ing them^®° showed that the these terms do not change 
the results qualitatively. The pairing interactions for the 
singlet and triplet channels are given by 

V'{q,ivi) = U + V{q) + ^U\s{q,^l) 

_ (i[/2 + 2UV{q) + 2V{qf)xc{q,uJi), (18) 

V\q,oJi) = V{q)-'^U\M,uJi) 
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Fig. 18. Phase diagram of the large- t-V model in tlie (T/t,V/t) 
plane sliowing competition between metallic, superconducting 
(SC), and charge ordered phases. The symmetry of the 
Cooper pairs in the superconducting phase is d^y , which is found 
near the quantum critical point {V/t)ci separating the metallic 
and charge ordered phases. [Reprinted figure with permission 
from J. Merino and R. H. McKenzie: Phys. Rev. Lett. 87 (2001) 
237002. Copyright (2001) by the American Physical Society.] 

_ ( 1 [72 + 2UV{q) + 2V{qf)xc{q, ^l), (19) 

respectively, where V{q) = 2t^(cos Qx + cos qy + cos{qx + 
Qy)) and uji is the Matsubara frequency (x and y- 
directions are taken along the p bonds in Fig. 12). Xs 
and Xc are the spin and charge susceptibilities, respec- 
tively. They are calculated within RPA as, Xs{q,^^i) = 

xoiq, Loi)/{i-Uxo{q, wi)), Xc{q, w;) = xoiq, u;i)/{i+{u + 

2V{q))xo{q,t^i)), where xo is the bare susceptibility. We 
note that the terms proportional to Xc in eqs. (18) and 
(19) represent effective pairing potentials due to charge 
fluctuation. This charge fluctuation contributes equally 
to V and to since it comes from the charge degrees 
of freedom. 

To determine the onset of the SC state, the linearized 
Eliashbcrg's equation in the weak coupling theory is 
solved. The obtained phase diagram on the {U, V) plane 
is shown in Fig. 19 at a fixed temperature T = 0.01. 
Aig{s*) represents the spin-singlet SC state with Aig 
symmetry, which is stabilized in the vicinity of both SDW 
(due to the nesting of Fermi surface) and CO (the three- 
fold state discussed in § 3.3) instabilities. When ^ ~ 0, 
the momentum dependence of the SC order parameter, 

A(fc), becomes dxy-li^e, similarly to the square lattice 

case.123.159.160 

Near the CO instability, on the other hand, a spin- 
triplet SC with symmetry is stabilized in addition to 
the spin-singlet state, although the eigenvalue is slightly 
smaller than that of the singlet pairing. The momen- 
tum dependence of the order parameter A(fe) shows 
that it is an /-wave pairing state since it changes the 
sign six times on the Fermi surface, which is similar to 
the isotropic triangular lattice case.^^^ This SC state is 
stable because the momentum dependence of V* gives 
large attractive interactions at the wave vectors such as 
Q = (27r/(3a),27r/(3a)). Actually Xc has a peak at Q 
when {U = -itpjV = 1.5tp) leading to a peak in V^{q). 
We note that, similarly to the square lattice case, this 
momentum Q is not due to the nesting instability of the 
noninteracting Fermi surface: Xo{q) has no notable peak. 




Fig. 19. A wealc coupling pliase diagram for the quarter-filled ex- 
tended Hubbard model for the 9-typc structure, on the (U, V) 
plane at T = 0.01 (tp is set as unity and tc = 0.4). i'^* The 
dashed lines correspond to the CO and SDW instabilities. The 
superconducting order parameter with Bgu symmetry is stabi- 
lized near the CO instability, although its eigenvalue is slightly 
smaller than that of the Aig symmetry. 



and Q is determined by the momentum dependence of 
V{q). A recent variational Monte-Carlo calculation^^® 
have also found the /-wave SC state to be stabilized, 
which suggests that the SC phase is near but not next 
to the CO phase. 

The stabilities of the dxy-wave pairing in the square 
lattice and the /-wave pairing in the anisotropic trian- 
gular lattice can also be understood from a naive real 
space picture. The nearest neighbor inter-site Coulomb 
repulsion repels electrons from the nearest-neighbor sites. 
The amplitude of the order parameter in real space be- 
comes larger at the four next-nearest-neighbor sites in 
the square lattice, while in the anisotropic triangular lat- 
tice there are six next-nearest-neighbor sites with large 
amplitude. These corresponds respectively to the dxy and 
the /-wave pairings. 

Let us briefly mention about SC found in Q;-(ET)2l3 
under uniaxial pressure which appears inside the CO 
phase; the SC phase transition takes place at a tempera- 
ture below an upturn of the resistivity is observed, which 
is probably due to CO.^^'' RPA calculations have been 
performed by Kobayashi et al.,^^'^ based on the recon- 
structed Fermi surface under the the horizontal stripe- 
type CO in the metallic phase obtained in the MF re- 
sults mentioned in 3.3. The results indicate that the SC 
phase can appear inside the middle CO metallic phase, 
in which small hole-pockets and electron-pockets both 
exist. ^''^ The suggested symmetry of this SC state is s- 
wave, of which the pairing interaction is mainly given by 
the spin fluctuation. This implies that in this case the 
charge fluctuation is suppressed due to the actual exis- 
tence of CO, and the pairing instability arises from the 
newly formed Fermi surface. 

4. Related Topics 

In this section we briefly mention several issues re- 
lated to CO in A2B systems that we have discussed. 
There are many transition metal oxides showing CO 
whose patterns can be identifled as the Wigner-crystal 
type one. Well known examples are the perovskite type 
compounds,^'^'^ such as 3D AMO3 or quasi-2D A2MO4, 
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with M being a transition metal. When A is half substi- 
tuted by another ion with different valence, as ^1/2^1/2' 
M can become mixed valent of half integer, namely, a 
"quarter" -filled situation is realized. In such cases, an 
"NaCl"-type CO for the 3D, or the checkerboard-type 
CO for the quasi-2D compounds is frequently observed. 
This CO state is actually behind the scene for the well- 
known coUosal magnetoresistance effect (CMR) in per- 
ovskite manganites.^^'' However, whereas the charge pat- 
terns reminds us of the CO states in the molecular sys- 
tems, the driving force in these systems is usually rather 
involved. Orbital degeneracy in the rf-electron under cu- 
bic/tetragonal crystal field, leading to the Jahn- Teller 
effect which is especially strong for the eg electron, fre- 
quently results in orbital ordering and makes the elec- 
tronic state highly anisotropic. This together with strong 
electron correlation give rise to interesting, but at the 
same time complicated phase diagrams. 

Another transition metal oxide compound has been 
discussed to show a similar CO transition as in the molec- 
ular conductors: NaV205, which has been intensively 
studied both theoretically as well as experimentally. 
The transition first found in the magnetic susceptibil- 
ity,^^^ from a behavior of paramagnetic ID localized spin 
systems to a spin-gapped state at 35 K, was revealed to 
be due to the CO transtion.^^° The average valence of 
the vanadium V^'^+ produces a quarter-filled dxy-ha.nd, 
which leads to an effective 2D quarter-filled model. It is 
the EHM, Hehm in eq. (1), on the so-called trestle lat- 
tice, where two-leg ladders along the a-axis are coupled in 
ilong the 6-axis. A MF study^^° predicted 
a CO state with a "zigzag" pattern, which is now con- 
firmed experimentally.^''^ As discussed in § 3, such a 2D 
model is difficult to study in a controlled way, and some 
authors studied the quarter-filled EHM on a two-leg lad- 
der system. ^^^'^''^ In fact, ladders are also found to be 
realized in molecular systems, ^''^ which is an interesting 
target for future studies. 

Recently, a triangular lattice system Naa;Co02 is at- 
tracting interest due to a SC state appearing when in- 
tercalated by water, H20.^'''' At a; = 0.5, a stripe-type 
CO state is stabilized, which is discussed to be coupled 
to the ordering of Na"*" ions.^''"'' This reminds us of the 
CO state in 6l-ET2RbZn(SCN)4 where the CO transition 
couples to the lattice degree of freedom and relaxes the 
geometrical frustration by lowering the symmetry of the 
lattice. Similar theoretical approaches to those discussed 
in this review have been applied to the EHM^''^ and the 
"t-V" model'^'"'' appropriate for this compound. 

In the research field of molecular A2B materials, one 
important direction of evolution is the inclusion of d- 
elcctron spin introduced in the B unit coupled to the 
organic 7r-electron on A, the so-called n-d systems. ^''^ 
The role of CO here is not clear as it seems that most 
of such materials are realized in effective half-filled sys- 
tems such as in the k- and A-type structures. It is desired 
that TT-d systems with a quarter-filled 7r-band would be 
investigated. It is to be noted that a theoretical work on 
slightly different molecular system TPP[FePcCN2]2 (Pc: 
pentacene), where the 5=1/2 Fe^+ ion is implanted 
in the donor tt molecule itself, showed that the interac- 



tion between localized spins would highly stabilize the 
CO state and predicted a novel ferromagnetic ground 
state. ^''9 

Many different kinds of molecular conductors are con- 
tinuously and constantly synthesized. We expect that CO 
states will appear ubiquitously in the new materials as 
well. Sometimes the effects of CO may be secondary; one 
recent example is the case of AB^ with x being an incom- 
mensurate value close to 1/2, i.e., with an off quarter- 
filled band. An "incommensurate Mott insulator" has 
been theoretically proposed, possibly realized in re- 
cently synthesized (MDT-TS)(Aul2)o.44i/*^ where the 
CO instability at quarter-filling indirectly controls the 
peculiar Mott transition of the system. 

5. Summary 

To summarize, we have reviewed theoretical stiidics 
on charge ordering and related phenomena in the 2:1 
charge transfer molecular conductors expressed as A2B. 
The charge ordered states are successfully described by 
extended Hubbard models which starts from the tight- 
binding model treating the anisotropy carefully enough 
and consider the electronic correlation, i.e., not only 
the on-site but also the inter-site Coulomb repulsion, 
and in some cases with the coupling to the lattice. 
The physics therein is rich, revealed by considerably nu- 
merous amount of studies, but still many issues arise 
from extensive experiments which anticipate future in- 
vestigations. Encouraged by this success of understand- 
ing the apparently complex system based on construct- 
ing microscopic models and treating them by different 
theoretical techniques, frontiers of research in this field 
are just about to expand. For example, recent activi- 
ties on the photo-induced phase transitions in A2B sys- 
tems^^'"'' requires theoretical efforts with new methods 
to explore their nonequilibrium dynamics from the mi- 
croscopic models. Furthermore, there is ambitious pro- 
posals^*"^ to apply our knowledge to even more complex 
molecular assemblies such as bio-related materials. 
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